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ABSTRACT 

Using relativistic Hartree-Fock approximation, we investigate the properties 
of the neutron-star matter in detail. In the present calculation, we consider 
^ not only the tensor coupling of vector mesons to octet baryons and the form 

CLf factors at interaction vertexes but also the internal (quark) structure change of 

baryons in dense matter. The relativistic Hartree-Fock calculations are performed 
in two ways: one is the calculation with the coupling constants determined by 
SU(6) (quark model) symmetry, the other is with the coupling constants based 
on SU(3) (flavor) symmetry. For the latter case, we use the latest Nijmegen 
(ESC08) model. Then, it is very remarkable that the particle composition of 
the core matter in SU(3) symmetry is completely different from that in SU(6) 
symmetry. In SU(6) symmetry, all octet baryons appear in the density region 
below ~ 1.2 fm~'^, while, in the ESC08 model, only the S~ hyperon is produced. 
Furthermore, the medium modification of the internal baryon structure hardens 
the equation of state for the core matter. Taking all these effects into account, we 
can obtain the maximum neutron-star mass which is consistent with the recently 
observed mass, 1.97 ± 0.04Mq (PSR J1614-2230). We therefore conclude that 
the extension from SU(6) symmetry to SU(3) symmetry in the meson-baryon 
couplings and the internal baryon-structure variation in matter certainly enhance 
the mass of neutron star. Furthermore, the effects of the form factor at vertex 
and the Fock contribution including the tensor coupling due to the vector mesons 
are indispensable to describe the core matter. In particular, the Fock term is very 
vital in reproducing the preferable value of symmetry energy, (~ 30—40 MeV), 
in nuclear matter. 
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1. Introduction 

Neutron stars are composed of the densest form of hadrons (baryons and mesons) and 
leptons. However, the detail of particle fractions and properties of the core of neutron star 
are not yet theoretically understood. Observations of the mass and/or radius of neutron star 
can provide the stringent constraint on the equation of state (EoS) of dense matter in the 
core. The best determined pulsar mass is 1.4414 ± O.OOO2M0 (the Hulse- Taylor pulsar), and 
masses of most other pulsars are found to be close to this canonical value. However, several 
neutron stars with heavier masses have recently been discovered. Radio timing observations 
of three Post Keplerian parameters led to the most preci se measurernent o f the mass of a 
milhsecond pulsar of 1.667 ± O.O2IM0 (PSR J1903-0327) JPreire et al.lboilh . Furthermore, 
Shapiro delay measurements from radio timing observations of the binary m illisecond pulsar 



(PSR J1614-2230) have indicated a mass of 1.97±O.O4M0 of the neutron star (jPemorest et al. 
2010h . 



Nuclear forces play a very important role in the core of neutron star, acting in concert 
with gravitational forces to form a compact object. The observed mass of a neutron star 
cannot be understood without counting the pressure arising from the strong nuclear force, 
which opposes gravitational collapse. The EoS for the core matter has been firstly inves- 
tigated in terms of only nucleons (A^) and leptons (£). Such a EoS can be very stiff, and 
thus yield the maximum allowable mass of a neutron star such as Mm nx — 2.8Mf;^, which i s 
close to the absolute upper bound on Mmax stemming from causality ( iHaensel et al.l 120071 ). 
However, in dense matter, hyperons {Y) generated through weak interaction may also par- 
ticipate in the EoS. This is true even at low density, in the region around twice the normal 
nuclear matter density, ?t.o(= 0.15fm~^), because some models predict the low threshold den- 
sity for the creation of hyperons in beta equilibrium. Therefore, there is no reason to deny 
the appearance of hyperons in dense matter. When hyperons are considered, the inevitable 
softening of the EoS occurs, which irnplies that the maximum neutron-star mass becomes 
Mmax ^ I.SMq ( Isee. e.g. Burgio et al.ll201ll : IVidana et al.ll201ll ). Such a low maximum mass 
is only marginally consistent with the Hulse- Taylor pulsar, but is already contradicted by 
the pulsars of PSR J1903-0327 and PSR J1614-2230. 

To explain the EoS for a neutron star, niany p eople have recently used Quantum Hadro- 



dynamics (QHD) (Isee. e.g. Serot fc Waleckalll984[ ) in relativistic mean- field (RMF) or rela- 



tivistic Hartree (RH) approximation, in which the baryons are treated as point-like objects 
and several non-linear, self-interaction terms of the scalar (a) and/or vector (w, p, 0) mesons 



a heavy neutron-star mass ( 


Todc 


-Rutel & Piekarewicz 


2005; 


Dexheimer & Schramm 


2008 


Shen et al. 


2011; 


Bednarek et al. 


2OI2I 


Weissenborn et al.l 


2012a 


)• 



- 3 - 



In general, an effective field theory at low energy will contain an infinite number of 
interaction terms, which incorporate the compositeness of hadrons, and it is then expected 
to invol ve numerous couplings wh ich may be nonrenormalizable. Naive- dimensional analysis 
(NDA) (IManohar &: Geor gill 19841 ) can provide an organizing principle to make sensible cal- 
culations for su ch system. For the nuclear many-body problem , the power counting scheme 
based on NDA (IFurnstahl et al.lll997l : iFurnstahl fc Serotll2000l ) allows us to expand the La- 
grangian (or the total energy) in terms of the ratios of scalar and vector densities to the 
factor f^A (the pion decay constant, and the large mass scale, A ~ IGeV). Such ra- 
tios are usually between 1/4 and 1/7 around uq, which can certainly serve as expansion 
parameters. From the point of view of the power counting, the tensor force due to the vr- 
and p-meson exchanges, which plays a very crucial role in nuclear forces, can, however, be 
viewed as a higher-order correction in a nuclear matter, and thus the tensor contribution to 
the total energy is not large in usual nuclei. However, in the core of neutron star where the 
density can grow up to ~ lOno, the expansion supported by the power counting is no longer 
valid, because the convergence becomes worse with increasing density. Thus, it is interesting 
to study how the tensor interaction affects the properties of the dense core matter. To take 
it into account, it is necessary to consider the Fock term, as well as the Hartree term. The 
Fock contribution also ensures the effect of anti-symmetrization for baryon wave functions 
in matter, which is not included in the RH calculations. Therefore, it is highly desirable to 
calculate the EoS within relativistic Hartree-Fock (RHF) approximation. 

At sufficiently high density, it seems unlikely that hadrons will maintain their identity as 
confined entities of quarks and gluons, but the interaction between overlapping hadrons will 
result in partial or full deconfinement. Even around no, it is well recognized that the value of 
the quark cond ensate, (qq), is not t he same as that in vac uum, and that it decreases by about 
20 — 30% at no ( ICohen et al.lll992t iTsushima et al.ll2007| ). Because the quark condensate has 
Lorentz- scalar character, the decrease in (qq) directly modifies the quark mass and thus leads 
to the change of quark wave function in a hadron. It implies that the hadron properties at 
finite density are certainly different from those in vacuum. Thus, even in the density region 
below the deconfinement of quarks and gluons, the EoS should contain all this physics. 



In the quark- meson coupling (QMC) model (iGuichonI Il988l : ISaito fc Thomaslll994al ). 
the properties of nuclear matter can be self-consistently calculated by the couphng of scalar 
(a) and vector (u) fields to the quarks within the nucleons, rather than to the nucleons 
themselves. As a result of the scalar coupling to the quark, the internal structure of the 
nucleon is modified with respect to the free case. In particular, the decrease of the quark 
mass caused by the attractive scalar field reduces the quark scalar density in the nucleon. 
Because the quark scalar density itself is the source of the a field, this provides a mechanism 
for the saturation of nuclear matter, where the quark structure plays a vital role. 
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In the past few decades, the QMC model has been e xtensively deve l oped and applied 



to various nuclea r phenomena with tremendous success ( iGuichon et al.l Il996l : ISaito et al. 



19961 . 119971 . 120071 ). In fact, the evidence for the medium modification of nucleon structure 



in a nucleus was observed in polarization transfer measurements in the quasi-elastic (e, e'p) 
reaction at th e Thornas Je fferson Laboratory, and the result supports the prediction of the 
QMC model (jBrooksll201ll ). It also seems vital to consider the internal str ucture change of 



the riucleon to understand the European Muon C ollaboration (EMC) effect ( iGeesaman et al. 
I995I : ISaito fc Tho"^ll994bt Icioet et allbooel l. 



The QMC model has recently been extended to include the quark-quark hyper fine in- 
terac tions due to the exchanges of gluon and pion based on chiral symmetry (iNagai et al. 
20081 ). and it has been applied to hyperons in a nuclear medium ( iMiyatsu fc Saitd l2009l ). 
We call this new version the chiral QMC (CQMC) model. The hyperfine interaction plays 
an important role in the in-medium baryon spectra. In particular, it can explain why the A 
feels more attractive force than the E or S in matter. The QMC model can offer a unique 
opportunity to self-consistently investigate the composition of the core of neutron star in 
terms of the quark and gluon degrees of freedom. 



In the previous calculation by lMiyatsu et al.l (|2012[ ). using the QHD, QMC and CQMC 
models, we have briefly reported the effects of Fock term, tensor coupling of vector raesons 
and baryon structure va riation on the EoS for the core matter (jsee also Stone et al.l 120071 : 
Whittenbury et al.ll2012l ). In the present paper, we fully investigate how the Fock term 
including the tensor interaction and the baryon structure change contribute to the EoS 
in RHF approximation. We also study the effect of f orm factors at inte raction vertexes. 
Furthermore, using the latest Nijmegen (ESC08) model (jRijken et al.ll2010l ). we examine the 
extension from SU(6) (quark model) symmetry to SU(3) (flavor) symmetry in the meson- 
baryon coupling constants. This extension leads to the significant changes in the particle 
fractions in the core matter and thus in the maximum mass. Our findings are summarized 
in the last section. 



2. Relativistic Hartree-Fock calculation 



We present the relativistic formulations of QHD and QMC models in Hartree-Fock 
approximation . 
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2.1. QHD with non-linear a terms 

In QHD, the octet baryons (proton (p), neutron (n), A, 2°^) and mesons (cr, u, 

TT, p ) are assumed to be structureless. Then, the Lagrangian density is given by 

£ = Cb + Ce + Cm + Ant, (1) 

where 

'Cb= Y1 ^B{il^.d^ - MB)iJB, '^i= Mh^.^>' - (2) 

with ipB^i) the baryon (lepton) field and Mb {fni) the baryon (lepton) mass in vacuum. The 
sum B runs over the octet baryons, and the sum i is for the leptons, e~ and fi~. For the free 
baryon masses, we take Mn = 939 MeV, Ma = 1116 MeV, Ms = 1193 MeV and Ms = 1313 
MeV, respectively. 

The meson term reads 

+ \mlp,-p^-^R,,-R^^ + ^{d,n-d^n-ml7r'), (3) 

with 

W^^ = df,uj^ - d^Uf,, Rf,y = d^p^ - d^p^, (4) 

where a and u are isoscalar mesons, while vf and p are isovector mesons. The meson masses 
are respectively chosen as mo- = 550 MeV, m^j = 783 MeV, = 138 MeV and nip = 770 
MeV. 



The interaction Lagrangian is given by 

Cint = y^v^ 



B 



fioB 



gpBi.p" ■ Ib + IS^Mv^V^' ■ Ib - —1,1,0^^ ■ Ib 

Zjvl m 



i^B, (5) 



where Ib is the isospin matrix for baryon B (we set Ib = when B is iso-singlet) and th e 
common, scale mass A4 (m) is taken to be the free nucleon (pion) mass (jRijken et al.ll2010l ). 
The cr-, u-, p-, n-B coupling constants are respectively denoted by g„B, QujB, QpB and /jts, 
while fujB and fpB are the isoscalar-tensor and isovector-tensor coupling constants. In RHF 
approximation, the meson field are replaced by the constant mean-field values: a, u and p 
(the p° field). Note that the mean-field value of the pion vanishes. 
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In the present calculation, we add non-linear (NL) terms of the a, 

Jr(-^ 92-3 , 93-4 (a\ 

U{cr) = ja + , (6) 

to the Lagrangian density, because the EoS given by the naive QHD is too hard, i.e., 
the incompressibility, K^, for symmetric nuclear matter is calculated to be ~ 550 MeV, 
which is apparently too large in comparison with the experimental value, 210 — 300 MeV 
( iDanielewicz et al.ll2010| ). Here, g2 and are additional coupling constants. We call this 
model QHD+NL. 

Although the baryons are assumed to be structureless in QHD, the effect of the size can 
be included if we introduce a form factor at each interaction vertex. In fact, at the Hartree 
level, it is unnecessary to include it because the momentum transfer in the meson exchange 
between two baryons vanishes. However, because the exchanged momentum can be finite in 
the Fock term, the effect of form factor may become significant as the density grows up. We 
choose a dipole function as the form factor: 

^^(^') = TT^TTi^' 

(1 - p^Af) 

where is the (four) momentum transfer, Aj is a cutoff parameter and i specifies the 
interaction vertex which will be discussed below (see also Table [T]). Note that the effect 
of form factor diminishes with increasing Aj. To study the effect of baryon size, in the 
interaction Lagrangian density, we replace all coupling constants (except for g2 and gs in 
Eq.(|6])) with those multiplied by the form factor, Eq. ([7]). 

To sum up all orders of the tadpole (Hartree) and exchange (Fock) diagrams in the 
baryon Green's function, Gb, we use the Dyson's equation 

GbW = G%{k) + Gl{k)^B{k)GB{k), (8) 

where fc^ is the four momentum of baryon, is the baryon self-energy and is the 
Green's function in free space. In a matter, the baryon self-energy is generally written as 

^B{k) = J:%{k) - 7oS°j(A;) + (7 ■ k)^Uk), (9) 

with k the unit vector along the (three) momentum k and the scalar part (the time 

component of the vector part) [the space component of the vector part] of the self-energy. 
Furthermore, the effective baryon mass, momentum and energy in a matter are respectively 
defined by 

M*(fc) = MB + J:%{k), (10) 
^b' = {k*B^k*^) = {k' + ^Uk),k + ki:i{k)), (11) 

EUk) = fe + M*2(d'^'. (12) 
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The baryon self-energies in Eq. (|9]) are then calculated by ( iBouyssy et al. I ll994i iMiyatsu et al.l 



20121) 



-gaBCr 



(I. 



BB' 



B',i ^ ' 
-QloB^^ - gpB{lB)3p - 



dq q 

{IhB'f 



2E*^,{qy 



B' 



7l (4vr)2A: J, 



dq qAi{k,q), 



i^BB') 

il Jo 



dqq 



'A:,C.{k,q) + ^^D,{k,q) 



2E*^,{q) 



(13) 
(14) 
(15) 



where kp^ is the Fermi momentum for baryon B and the factor, /^^/, is the isospin weight 
at the meson-i?i?' vertex in the Fock diagram. In the present calculation, we ignore the 
retardation effect (i.e., the energy transfer between two inte racting baryons vanishesjll, whic h 
gives at most a few per cent contribution to the self-energy (Isee. e.g. Serot fc Waleckalll984j ). 
The index i in the sum and the functions Ai, Bi, Ci and Di in Eqs. (fT3|) - (fT5|) are explicitly 
specified in Table [H in which the following functions are usecj^: 



e,(fc,g) 

^i{k,q) 

^i{k,q) 
U,ik,q) 
^i{k,q) 



A? 



1 



In 



M+{k,q)L-{k,q) 



Akq 



M-ik,q)Ltik,q) 
[{k' + q' + m^) e,{k,q)-AfNf{k,q 



71=1 



{e + q^- m^/2) <!>i{k, q) - kqOi^k, q) + ^(fc, g), 

+ g2) q) _ kqQi{k, q) + ^(A;, g), 

[fce,(/c,g)-2g$,(A;,g)], 



where 



Vti{k,q) 

Lf{k,q) 
Mt{k,q) 

Nnk,q) 



A? 



[NUk,q) + ie + q' + A^)Nf{k,q 



Akq 
K} + {k±qf., 

m} + {k± qf 



n 



[Lt{Kq)r-[L.{k,q)]- 



(17) 

(18) 
(19) 
(20) 



(21) 

(22) 
(23) 

(24) 



^In this approximation, the form factor is given as Fi{p^) ^1 +p^/Af j , where p is the (three) 
momentum transfer. 



^In the loop integral, we have removed the so-called contact interaction term (jBouyssv et all Il994 : 



Krein et al.lllQQOl) . 
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The mean-field values of u and p are respectively given by the usual forms 

y'^PB, p=y'^(4)3PB, (25) 



B 



B 



where the density of baryon B is ps = C^Q^t^ )^FB wi^h Jb the spin of B. 

On the ot her hand, combining wi th Eqs. (IT3l)- (fT5l) . the a value should be self-consistently 
calculated by dWeber fc Weigellll989l ) 



a 



E9ctB s 1 / -2 , -3 
—Pb 2y92(^ +93(^ 
g mi 

where the scalar density of baryon B reads 



(26) 



Pb 



2Jb + 1 
27r2 



dk k 



Em 



{21) 



The total energy density and pressure can be divided into the hadron and lepton parts, 
namely e = + and P = Ph + Pe, where the hadron energy density, en, and pressure, 
Ph, include both the baryon and meson contributions. Then, they are expressed as 



B 



2Jb + 1 

JO 



dk 



TB{k) + \vB{k) 



(92 . 93- 

— cr 

2 V3 2 



(2^ 



with 



TB{k) 



MBM*g{k) + kkl 

Em ' 



^ M%{k)T.%{k) + k*j,i:m iy, 
^b[K) - ^Blfcj, 



E*m 



and 



d 



^B 



(29) 



(30) 



where the total baryon density, n^, is given hy ub = J2bPb- These energy density and 
pressure give the EoS for the core matter of neutron star. 



2.2. QMC and CQMC models 

In QHD, QcrB is constant at any density because the baryons are point-like objects. In 
contrast, in the QMC and CQMC models, the a-B coupling constant, (^o-s, has the a- field 
dependence, which refiects the baryon structure variation due to the interaction in matter. 
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However, the mediu m effect does not af fect tlie coupling of vector meson to baryon in mean- 



field approximation (ISaito et al.l 120071 ). Furthermore, in a uniformly distributed matter, 
it may not be necessary to consider the internal structure variation of mesons explicitly, 
because the meson field is no longer dynamical hnt just auxiliary, namely the meson fie ld (or 



the m ean-field) can be replaced by the baryon density or the baryon scalar density (jSaito 



2OIOI). 



In the QMC model, the MIT bag model is simply used to describe the baryons, in which 
three quarks move freely. As an example, we consider the nucleon mass in matter. When we 
calculate the mass in RH approximation, it is given as a function of the scalar mean-field 

M^(a) = M;v-(7<xiv(a)a, (31) 

because the scalar field couples to a quark in the nucleon and modifies the quark mass in 
matter. Here, gaN^o') is the a-N coupling constant, which depends on the a field, and will 
be discussed further below. 

Because the value of a is small at low density, the effective nucleon mass in matter, M^, 
can be expanded in terms of a as 

In the QMC model, the interaction Hamiltonian between the nucleon and the a field at the 
quark level is given by Hint = —^9% Jy df ipgaipq, where is the cr-quark coupling constant, 
ipq is the quark wave function and V is the nucleon volume, and the derivative of with 
respect to a is given by 

(^) = X "^"^ ^ -USr,{a). (33) 

Here we have defined the quark-scalar density in the nucleon, S'Ar(a), which is itself a function 
of the scalar field, by Eq.( l33|) . Because of a negative value of f-g^j, the nucleon mass 
decreases in matter at low density. 

Moreover, we define the scalar-density ratio (or the scalar polalizability), SNicr) / Sn{0), 
to be Cj\f{a) and the a-N coupling constant at a = to be g^-N (i-e., g^N = gaNip- = 0)): 

C^{a) = 5^(a)/S;v(0) and g,r^ = 3(?^S;v(0). (34) 

Comparing with Eq.f l3ip . we find that 

_ _g^^CN{^) = [gaN{o-)a] , (35) 



oa J da 
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and that the nucleon mass is 



Mn - gaN<y - 25'<TivC^(0)(T^ + 



(36) 



In general, Ctv is a decreasing function because the quark in matter becomes more relativistic 
than in free space. Thus, C^(0) takes a negative value. If the nucleon were structureless, Cn 
would not depend on the scalar field, that is, Cjq would be constant {C^ = 1). Therefore, 
only the first two terms in the right hand side of Eq. fl36|) remain, which is exactly the same 
as the equation for the effective nucleon m ass in QHD. T his decrease in Cat constitutes a new 
saturation mechanism for nuclear matter (jGuichonlll988l ) , and very much reduces the nuclear 
incompressibility K^, compared with the value of QHD. Therefore, in the QMC model, it is 
not necessary to include the non-linear terms of a. 

In actual numerical calculations, we found that t he scalar-density ratio, CAf(cr), decrease s 
linearly (to a very good approximation) with gaNO' ( iGuichon et al.lll996l ; ISaito et al.l 119971 ). 
Then, it is very useful to have a simple parameterization for Cn: 



C^(cr) = 1 - Oat X (g^Ncr). 



(37) 



As a practical matter, it is easy to solve Eq. fl35l) for gaN{(^) in the case where C{a) is hnear 
in gaNO', as in Eq.( l37|) . Then one finds 



-{gaN<y] 



(38) 



M"^ = Mn - gaNicr)a = Mn - g^N 

so that the effective a-N coupling constant, g^Ni^), decreases at half the rate of CAr(a). 

For the octet baryons, this linear approximation is also applicable to the scalar-density 
ratio and the coupling constant in a nuclear medium: 



Cij(cr) = 1 - as X {g^NCr 
gasio) = g^B 



1 I _^ 

1 - Yl^aJVfT) 



(39) 
(40) 



with Ob a parameter and g^B the a-B coupling constant at ub = 0. This idea leads t o a 



new, simple scaling relationship among the in- medium baryon masses (jSaito et al.l 119971 ). 



Two quarks in a baryon can interact through the one-gluon exchange (OGE) (IGuichon et al 



20081 ) . In addition to it, the pion-cloud effect, which is governed by chiral symmetry and its 
small breaking, may also affect t he properties of baryons in a nuclear medi um. The CQMC 
model involves both the effects (iNagai et al.ll2008l ; iMiyatsu &: Saitd l2009l ). The hyperfine 
interaction due to the OGE and pion cloud is very significant in the in-medium baryon spec- 
tra. In particular, at the quark mean-field level, it splits the masses of A and S in matter, 
and, as it should be, the A feels more attractive force than the S or H does. 
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In the case of the CQMC model, we need two parameters, as and bs, to describe the 
(T-B couphng constant: 



gaB[(y) = QuBbB 



(41) 



while the scalar-density ratio can again be given by Eq. f l39p . This linear approximation is 
quite accurate up to ~ lOrio- The values of as and 6^ are listed in Table O Using Eq. lHTl) . 
the octet baryon mass in RH approximation is then given by = Mb — gcrB{o')o', which 
is shown in Fig. [H If we set = and 6^ = 1, gcrB{(^) is identical to the a-B coupling 
constant in QHD. We should note that the a- field dependence in the coupling constant in 
Eg. (HTD can also be understood as th e effect of many-body correlations in nuclear matter 
JCuichon fc Thomasl[20oi : ISaitolboiol ) 

For the RHF calculation in the QMC or CQMC model, we use the formulation developed 
in section \2A[ in which, instead of the constant g^B, the field-dependent coupling, QaBi^-), 
in Eq. (HTj) is used at each a-B vertex. 



3. Numerical results 



3.1. SU(6) symmetry 

We focus on how the Fock term including the tensor coupling and the form factor 
at interaction vertex contribute to the EoS for a neutron star. There is large ambiguity 
in determining the meson- hyperon coupling constants. Then, the SU(6) (quark model) 



symm etry is often assumed to fix the couphng constants (IReuber et al. 
1998h : 



1994 : iHuber et al. 



9aN 
9pN 
fnN 



3 3 „ 

-guA = -guT: = iga-E, 
1 1 

3 
5 



-guiN — -gpT. — gpE, 



guiN 



gpA = 0, 



:guiA 



zguT. — 'igu 



^pN 



fuiN 
guiN 
fuT. 

guB 

fpN 

gpN 



: -0.12, 

1 + 2k^n = 

— 2 — K^TV = 

3.7, K,p\ 



/ttA = 0, 
fwA 



0.76, 



gu>A 



-1.88, 

gpA 



(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
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«:ps = ^ = -| + |«:p^ = 0.88, (49) 

«:pH = ^ = ----/€p^ = -1.94, (50) 
5ps 5 5 

where we take f^j^/An = 0.08, and k, is defined by the ratio of the tensor (/) to vector (g) 
couphng constants. 

To see the effects of the Fock term and form factor, we consider the following cases: 

(1) RH calculation in QHD-NL, 

(2) RHF calculation (no tensor coupling, no form factor) in QHD-NL, 

(3) RHF calculation with tensor coupling (no form factor) in QHD-NL, 

(4) RHF calculation with tensor coupling and form factor (full calculation) in QHD-NL, 

(5) full RHF calculation in QMC, 

(6) full RHF calculation in CQMC, 

(7) full RHF calculation with hyperon-potential fit in QHD-NL. 

For simplicity, in cases of (4) - (7), all cut off parameters in the form fa ctors are assumed to 



be the common value, = 0.71 GeV^ (IFrauenfelder fc HenleylllQQll ). Once go-N is fixed 



the a-Y couplings, g^y, are determined through the relation Eq. fH2l) . In cases of (1) - (6), we 
use such values for gcrv- On the other hand, the recent analyses of hypernuclei and hyperon 
production reactions infer that the A, E or S, respectively, feels the potential, Ua,t., ; " 



—30, +30, —15 MeV, in a nuclear medium (jSchaffner &: Mishustinlll996l : llshizuka et al.ll2008l ). 



Therefore, in case (7), instead of Eg. (1121) . we determine g^A, QaY, and (^o-s so as to fit the 
suggested potential depths at no. Apparently, this leads to the breaking of SU(6) symmetry 
in the coupling constants. 

In Table [31 we list the coupling constants and several properties of symmetric nuclear 
matter at hq. We have four, free parameters in QHD-NL: ^f^-Ar, g^j^, g2 and 5^3. The first 
two parameters, g^N and g^iN, are adjusted so as to reproduce the nuclear saturation energy 
(—15.7 MeV) at uq in symmetric nuclear matter. The last two constants, g2 and (73, are 
determined so as to generate the incompressibility K„ (=253 MeV) and the effective nucleon 
mass {M^/M^ = 0.79) calculated by the QMC model (case (5)). For case (7), the coupling 
constants except for g^y are identical to those in case (4), and we use g'^j>^/4:7[ = 1.97, 
glj^/4:7i = 0.38 and S'^h/^tt = 0.52, which produce the required hyperon-potential depths. In 
the QMC and CQMC models (cases (5) and (6)), we have only two parameters, gaN and 
OujN, which are again determined so as to satisfy the nuclear saturation condition. We note 
that the values of K^, and calculated by the QMC are different from those by the 
CQMC, which is caused by the hyperfine interaction due to the gluon and pion exchanges 
between quarks, as discussed in section 12.21 
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Comparing cases (1) - (3), we can see that, because the Fock terms, especially the terms 
involving the tensor couplings, contribute as an attractive force, the coupling constants, 
g^M and g^]^, are reduced in cases (2) and (3) to satisfy the saturation condition for nuclear 
matter. In case (3), the Fock contribution at Uq reaches about 35 (25)% of the scalar (vector) 
self-energy. In contrast, the inclusion of form factors reduces the Fock contribution, which 
is shown in the contents of the self-energies in case (4), and g^N and g^^N are accordingly 
enhanced. Because of the quark and gluon degrees of freedom, the scalar and vector self- 
energies in the CQMC model (case (6)) are more enhanced (deeper) than in QHD-NL, which 
generates the lower effective nucleon mass in matter and the larger nuclear incompressibility. 

We should notice that the Fock term plays a very important role in the symmetry 
energy, 04. The value of in case (1) is very small, while, the Fock term enhances 04 very 
much (see cases (2) - (6)). We present the density dependence of in Fig. [21 In the RH 
calculation (case (1)), it is necessary to enhance the p-N vector coupling by about 45% by 
hand to obtain the observed symmetry energy, 04 ~ 30 MeV. 

In Figs. [3] and m we respectively show the equation of state and the particle fractions, Yi 
{i = n,p, ■ ■ ■), in the core matter of neutron star, which are calculated under the conditions 
of the charge neutrality and the /3-equilibrium in weak interaction. As expected, the hard 
EoS can be given by the QMC or CQMC model, while the softest one is generated by case 
(2). In Fig. m we can see all octet baryons in the RH result (case (1)). The Fock terms 
(without the tensor couplings) push upward the threshold densities for hyperons (except 
for E~ and E°), see case (2). Furthermore, the inclusion of the tensor coupling (case (3)) 
makes the threshold densities of A and S° very high (actually, S° disappear in the figure), 
while that of the S~ is somewhat pulled downward. In case (4), we can see that such tensor 
effect is weakened by the inclusion of form factors. In case (5), comparing with case (4), 
the thresholds of hyperons are again pushed upward, which is caused by the effect of baryon 
internal structure variation in matter. The difference between cases (5) and (6) is generated 
by the hyperfine interaction between quarks in octet baryons. In case (7), as expected, the 
E hyperons do not appear at all below ub = 1.2 fm~^, which is caused by the adjustment of 
the coupling constant (^o-s to fit the repulsive potential depth. Us — +30 MeV. 

We respectively present the neutron-star mass as a function of the radius and the prop- 
erties of neutron star in Fig. |5] and Table HI To calculate the neutron-star mass, we solve th e 



Tolman-Oppenheimer-Volkoff (TOV) equation (jTolmanlll934j : IOppenheimer fc Volkovlll939l ). 
As expected from Fig. |3l the CQMC model gives the heaviest mass. In contrast, although 
the mass is low, the case (2) gives the neutron star whose central density is very high, namely 
it is very compact. In SU(6) symmetry, the maximum neutron-star mass gradually increases 
from case (2) to (4). 
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3.2. SU(3) symmetry 

To study the role of hyperons on the properties of neutron star, it is very important 



to ext end SU(6) symmetry to the more general SU(3) flavor symmetry (jWeissenborn et al. 



2012bl ). Restricting our interest to three quark flavors, up, down and strange, SU(3) symme- 



try in flavor space can be regarded as a symmetry group of strong interaction. Because we 



consider only the octet representation for the baryons {J^ = octet), the SU(3) invariant 
Lagrangian can be constructed using the matrix representations for the baryons, B, and 
meson nonet (singlet state. Mi, and octet state, Ms). The interaction Lagrangian is given 
as a sum of three terms, namely one stemming from the coupling of the meson singlet to the 
baryon octet {S term) and the other two terms from the interaction of the meson octet and 
the baryons - one being the antisymmetric [F) term and the other being the symmetric (D) 
term: 

Cnt = -98^2 [«Tr([5, Ms]B) + (1 - «)Tr({B, Ms}B)] - (7i4^Tr(55)Tr(Mi), (51) 

where gi and gs are respectively the coupling constants for the meson singlet and octet states, 
and a (0 < a < 1) denotes the F/{F + D) ratio, namely the ratio of the antisymmetric 
contribution, F, to the sum of antisymmetric and symmetric, D, ones. In SU(3) symmetry, 
all possible combinations of couplings are then described in terms of four parameters; the 
singlet and octet coupling constants, gi and gs, a and the mixing angle, 6, relating the 
physical isoscalar mesons to the pure singlet and octet counterparts. 

The extended-soft-core (ESC08) model by the Nijmegen group m ay at present be th e 



most complete meson-exchange model for the baryon-baryon interaction (iRijken et al.ll2010l ). 
It describes not only the N-N but also the Y-N and Y-Y interactions in terms of the meson 
exchange based on SU(3) flavor symmetry. In this subsection, we calculate the EoS using 
the meson-baryon coupling constants provided by the ESC08 model. 

In this model, for example, the couplings of w-meson to baryons are related through 

QujA _ 1 gu;E _ 1 - V3z tan 9^ 

g^N l + v32;tan6'„ g^iN l + v3ztan6't, 

where z = gs/gi- The isoscalar u and (p meson are described in terms of the pure singlet, 
|1), and octet, |8), states as 

\u) = cos ^„|1) + sin ^„|8), (53) 
10) = -sin^^ll) + cose^|8), (54) 

where the mixing angle is chosen to be 6^ = 37.50°, which is very close to the ideal mix- 
ing, tan^;'^^'*' = 1/^2, namely ^j,'^""^ = 35.26°. The ratio z is flxed to be 0.1949 in the 
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ESC08 model, while z = l/y/G = 0.4082 in SU(6) symmetry. Thus, the ratio z is much 
smaller than the SU(6) va lue, which helps to enhance the maximum mass of neutron star 
( Weissenborn et al. 2012bl ). We show the coupling constants (except for g^B and g^^N) and 



the cutoff parameters in the form factors in Table [5l It should be noted that, in the present 
calculation, we use the dipole function for the form factor, Eq.([7]), while, in the ESC08 
model, it is expressed by a gaussian function (for details, see the footnotes of Table [5]). 

Using these coupling constants and cutoff parameters, we here examine the following 
four cases: 

(a) full RH calculation in QHD-NL, 

(b) full RHF calculation in QHD-NL, 

(c) full RHF calculation in QMC, 

(d) full RHF calculation in CQMC. 

The properties of symmetric nuclear matter at no and the rest of coupling constants are 
shown in Tables [6] and [71 As in the case of SU(6) symmetry, the coupling constants, g^N 
and gujN, are determined by the nuclear saturation condition. In QHD-NL (cases (a) and 
(b)), the coupling constants for the non-linear terms, g2 and g^, are again determined so as 
to fit the incompressibility and the effective nucleon mass to the values [K^ = 249 MeV and 
M^/Mj^ = 0.78) given by the QMC model (case (c)). As discussed in subsection 13. H the 
a-Y couphng constant (see column (A) in Table [7]) is fixed so as to reproduce the observed 
hyperon-potential depth, Uy- 

We show the symmetry energy, 04, as a function of in Fig. El It should be again 
emphasized that in case (a) the symmetry energy is clearly insufficient because of the absence 
of the Fock term. In contrast, 04 is very enhanced in cases (b) - (d). In Figs. [TlandlHl we 
respectively present the equation of state and the particle fractions in the core matter. The 
structure of particle fractions is very simple and impressive in all cases, namely only the 
H~ appears below ub = 1.2 fm~^. In the EoS, we can see the bend of pressure line around 
6 ~ 550 MeV/fm'^, which corresponds to the threshold point of the S~ production. Beyond 
the bend, the pressure in the RH calculation (case (a)) is enhanced compared with the other 
cases (b) - (d). It may be caused by the fact that, in comparison with cases (b) - (d), the 
production of E~ in case (a) is relatively suppressed for high densities, and that the fraction 
of neutron is thus not reduced much. This is also consistent with the fact that the threshold 
density for the production in case (a) is relatively high (> 0.6 fm~^), comparing with 
the other cases. As seen in subsection 13.11 because the tensor coupling of vector meson pulls 
downward the threshold density of the S~, the reduction of the threshold density in cases 
(b) - (d) may be originated by the effect of tensor coupling. 

In Fig. [9l we present the neutron star mass as a function of the star radius. Because 
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of the hard EoS, the RH calculation gives the heaviest maximum mass around 2.0Mq. The 
precise results for the neutron-star properties are presented in column (A) of Table [HI The 
CQMC model gives the mass of I.93M08 which is also consistent with the recently observed 
mass, 1.97± O.O4M0 (PSR J1614-2230). Comparing the results of cases (b) and (d), we 
conclude that the change of internal structure of baryons in matter can generate the mass 
fraction of ~ O.IMq. Furthermore, because the difference between the QMC and CQMC 
models is caused by the hyperfine interaction between quarks in a baryon, such interaction 
is also important and produces the mass fraction of ~ O.OTMq. 

We here notice that, because the coupling constants, g^B and q^jn, are determined so 
as to satisfy the desired conditions, namely the saturation condition for symmetric nuclear 
matter and the observed hyperon-potential depths, SU(3) flavor symmetry in the coupling 
constants is now partly broken. However, the violation of SU(3) symmetry may naturally be 
expected, because the core matter of neutron star consists of octet baryons whose particle 
fractions are not constant but vary depending on the density. Furthermore, the effective 
coupling constants in matter may, in general, vary depending on the surrounding matter, if we 
treat the matter using mor e sophisticated models such as the (relativistic) Dirac-Bru eckner- 



Hartree-Fock calculations ( Isee. e.g. Krastev fc Sammarrucall2006l : iDalen et al.ll2007l ). etc. 



However, it may be interesting to examine the neutron-star mass again, restoring SU(3) 
symmetry. For example, we try to recalculate the neutron-star mass by reducing the cu-Y 
vector coupling constants by the factor x{i) = qIjn/q^n^ where (7^^ and g^j^ are, respectively, 
the oj-N coupling constants for case {i = a, b, c or d) (see Table and that in free space 
(i.e., the value given in the ESC08 model). We choose the tensor coupling, /^y, so as to 
keep K^y given in the ESC08 model (see Table [5]). This reduction may partly restore SU(3) 
symmetry and bring the present calculation close to that with SU(3) symmetry. We then 
flnd that the factors are as follows: x(a) = 0.689, x{h) = 0.641, x(c) = 0.626, x{d) = 0.683. 
Correspondingly, the a-Y coupling constant should be readjusted so as to reproduce the 
correct potential depth, Uy, which is shown in column (B) in Table [3 This recalculation 
gives the maximum neutron-star mass presented in column (B) of Table [HI namely the mass 
is reduced by about 8% comparing with the result in column (A), and thus the mass cannot 
reach the value of 1.97 ± 0.04Mq. Fig. [TU] shows the particle fractions, where not only the 
S~ but also the A hyperon appear in the core matter. The threshold density for the A 
production is ~ 0.4 — 0.5 fm~'^, which is very close to that for the H^. As in cases of column 
(A), their threshold densities in case (a) are somewhat higher than those in cases (b) - (d). 



•^In the previous calculation bv iMivatsu et al.l (j2012l ). we have reported Mmax — 2.02Mq in the CQMC 
model, in which the CQMC predicts Ky — 300 MeV because the form factor is not included. The difference 
between the previous and present results is mainly caused by the effect of form factors at interaction vertexes. 
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Therefore, the RH result again gives the heaviest neutron-star mass. 

4. Summciry and conclusions 

Using relativistic Hartree-Fock approximation, we have studied the properties of the 
core of neutron star in detail. In the present calculation, we have considered not only the 
tensor coupling of vector mesons to octet baryons and the form factors at vertexes but also 
the internal (quark) structure change of baryons in matter. 

We have performed the RHF calculations in two ways: one is the calculation with the 
coupling constants determined by SU(6) (quark model) symmetry, the other is with the 
couphng constants based on SU(3) (flavor) symmetry. Then, from the result of the former 
calculation, we have found the following important features. The Fock term (except for the 
tensor coupling) makes the EoS for neutron star rather soft, while the tensor coupling pushes 
the threshold of hyperon production upward and thus the EoS becomes somewhat hard. 
Because the form factor reduces the effect of exchange term and the scalar self-energies for 
baryons due to the tensor coupling, the a and u mean fields arc correspondingly enhanced to 
keep the reasonable values of Ky and the nucleon mass in matter. Then, the enhancement of 
u mean field hardens the EoS. Furthermore, the medium modification of the internal baryon 
structure also makes the EoS hard, and enhances the maximum neutron-star mass. If the a- 
hyperon coupling constants are readjusted so as to reproduce the observed hyperon-potential 
depths (although such change eventually breaks SU(6) symmetry), the particle composition 
in the core matter is considerably modified, namely the E hyperons disappear below ub = 1-2 
fm~^. However, the calculated maximum neutron-star mass in SU(6) symmetry cannot reach 
the massive neutron-star mass, 1.97± O.O4M0 (PSR J1614-2230). 

In the calculation based on SU(3) symmetry, we have used the coupling constants pro- 
vided by the latest Nijmegen (ESC08) model. Then, it is very remarkable that the structure 
of particle composition in matter is completely different from that in the case of SU(6) sym- 
metry, namely all hyperons except the disappear below Ub = 1-2 fm~^. Accordingly, the 
EoS becomes hard, and thus, in the RH or CQMC calculation, we can obtain the mass which 
is consistent with the observed mass of 1.97 ± O.O4M0. Because, in the present calculation, 
the coupling constants for the non-linear cr terms in QHD are determined so as to reproduce 
the values of and M'^ given by the QMC model, the coupling of the u meson to the nu- 
cleon in the RH case is rather strong compared with those in the RHF cases. Thus, the RH 
calculation gives the hardest EoS, which generates the heaviest neutron-star mass, 2.OM0. 
The medium modification of the internal baryon structure also plays the significant role in 
determining the neutron-star mass. Such effect generates the mass fraction of about O.IMq. 



- 18 - 



Because the coupling constants, gas and Quin, are determined so as to satisfy the nu- 
clear saturation and the hyperon-potential depths at uq, SU(3) symmetry in the coupling 
constants is partly broken in the present calculation. However, such violation may be jus- 
tified, because the ground (vacuum) state of the core matter consists of the octet baryons, 
which fill (positive) baryon-energy levels to the Fermi momenta, and each baryon fraction 
inside the core varies depending on the density. It implies that SU(3) symmetry is explicitly 
broken in the ground state. However, we have checked how the neutron-star mass is modified 
when SU(3) symmetry is restored. In fact, by reducing g^y by hand, which may (partly) 
restore SU(3) symmetry, we have recalculated the neutron-star mass. Then, we have found 
that the maximum mass cannot reach the desired value, 1.97±0.04Mq. Thus, how we deter- 
mine the coupling constants, especially the coupling constants of vector mesons to baryons, 
is very important to obtain a heavy neutron-star mass. 

In conclusions, we have found that the extension from SU(6) symmetry to SU(3) sym- 
metry in the meson-baryon couplings and the internal baryon-structure variation in matter 
certainly enhance the maximum neutron-star mass. Furthermore, the effects of the form 
factor at vertex and the Fock term including the tensor coupling due to the vector meson are 
indispensable to describe the core matter of neutron star. In particular, the Fock contribu- 
tion is very vital to reproduce the preferable value of symmetry energy, 04 (~ 30 — 40 MeV), 
in nuclear matter. Including all these effects, it becomes possible to explain the mass of 
PSR J1614-2230 pulsar. In the usual RH calculations with the coupling constants in SU(6) 



as well as the a meson, to obtain a heavy neutron-star mass ( 


Todd-Rutel & Piekarewicd 


2OO5I: 


Dexheimer & Schramm 


2008; 


Shen et allboil: 


Bednarek et al.ll2012l:lweissenborn et al. 


2012a) 



However, it should be emphasized that such self-interaction terms are not indispensable in 
the present RHF calculation. 

Finally, we give several comments on the future work. It may be important to consider 
the K or K* meson exchange between two baryons, because they generate the mixing of octet 
baryons in matter. For example, the mixing of A^-A, A-S, A^-S or S-A occurs through the 
K or K* meson exchange. Furthermore, the tt or p meson mixes the S-A hyperons. These 
mixing effects may change the baryon composition i n the core matter. It is also importan t 
to introduce temperature into the RHF formulation (ISaito et al.lll989l : ISoutome et al.lllQQOl ). 



In the present study, we do not consider the short range (two-baryon) correlation in 
matter, which may affe ct the EoS at high density. In fact, using the latest ESC08 model, 
Schulze fc RijkenI (120111 ) have already performed the Brueckner-Hartree-Fock calculation for 
the neutron- s tar ra a tter, and reported ve r y low maximurn mass es below 1.4M0 (see also, 
Akmal et al.l (Il998[ ): iNishizaki et al.l ( 120021 ); iLogoteta et al.l ( 12012| )). Because relativity must 
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be very crucial to investigate the high density matter, it is very desired to perform the 
(relativistic) Dirac-Brueckner-Hartree-Fock calculation including hyperons. 

At very high density, the quark and gluon degrees of freedom, rather than the hadron 
degrees of freedom, may take place in the core matter. Because the degrees of freedom in a 
quark- gluon matter is generally large, it is necessary to assume a r ather strong correlati on be- 



tween quarks and gluons to support a massive neutron-star mass (jBratovic et al.ll2012l ). Fur- 



thermore, in the crossover between the hadron and quark-gluon pha ses there may exist rich 



non-p erturbative structure such as color superconducting phases etc (IFukushima fc Hatsuda 



20111 ). It would be very intriguing to investigate how such degrees of freedom contribute to 



the EoS and the maximum mass of neutron star. 



The authors (especially, T. M.) would like to thank Myung-Ki Cheoun, Tomoyuki 
Maruyama, T. Kajino and T. Takatsuka for fruitful discussions on the EoS for neutron 
stars. 
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Fig. 1. — The CQMC result for octet baryon masses (as functions of the scalar field g^NO') 
in symmetric nuclear matter. We note that the difference between the results by the exact 
calculation and the parameterization, Eq. (HT]) . is indistinguishable. Here, we use the cou- 
pling constants, g^B, given by SU(6) symmetry, namely g^N = fs'^A = ffi'^s = ^g^rE (see 
section 13. ip . 
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Fig. 2. — Symmetry energy 04 in SU(6) symmetry. 
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Fig. 5. — Neutron star mass versus radius in SU(6) symmetry. 
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Fig. 6. — Symmetry energy, 04, in the ESC08 model. 
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Fig. 7. — Equation of state (pressure versus energy density) in the ESC08 model. 
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Fig. 8. — Particle fractions, Yi, for cases (a) - (d). 
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Fig. 9. — Neutron star mass versus radius for cases (a) - (d). The shaded area shows the 
mass of pulsar (PSR J1614-2230), 1.97 ± O.O4M0. 
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Fig. 10. — Particle fractions, 1^, for cases (a) - (d). The uj-Y coupling constants are reduced 
to restore the SU(3) relation (for details, see the text). 
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Table 1: Functions Ai, B^, Ci, and Di. The index i is specified in the left column, where 
V{T) stands for the vector (tensor) coupling at each meson-BB' vertex. The last row is for 
the (pseudovector) pion contribution. 



(7 

UJvv 

ujtt 

OJvT 

Pvv 
Ptt 

PVT 



glB{(y)Qa 



\2M 



Bi 



-3 ih^Tml^u 



2 ( ^ ) n 



1 9 / AiB9uiS.'\ p 
-^^ I 2M ) ^' 



1 9 / fpBOpB \ p 

-^^ ' 2M r- P 



Table 2: Values of and 6^ for octet baryons in the QMC or CQMC model. 





QMC 




CQMC 






as (fm) 


bB 


qb (fm) 


bB 


N 


0.179 


1.00 


0.118 


1.04 


A 


0.172 


1.00 


0.122 


1.09 


E 


0.177 


1.00 


0.184 


1.02 




0.166 


1.00 


0.181 


1.15 
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Table 3: Coupling constants (in SU(6) symmetry) and calculated properties of symmetric 
nuclear matter at uq. The values of the incompressibility, K^, the symmetry energy, 04, and 
the self-energies for the nucleon, and S^, are in MeV, and, in the parenthesis of the 
self-energies, the left (right) number is for the Hartree (Fock) contribution. We note that 
the space component of the vector self-energy, S^, is very small (at most ~ 5 MeV). Note 
that, for case (7), the coupling constants except for g^y are identical to those in case (4). 



case 


9 a N 

47r 


9u,N 


92 (fm-i) 


93 




04 


Mis, 


ys 


yO 


(1) 


5.73 


5.66 


21.5 


34.1 


253 


21.7 


0.79 


-193 (-193,0) 


-134 (-134,0) 


(2) 


4.67 


3.91 


10.7 


170 


253 


33.3 


0.79 


-193 (-157, -36) 


-135 (-92,-43) 


(3) 


3.78 


4.23 


10.5 


235 


253 


35.4 


0.79 


-193 (-126,-67) 


-133 (-100,-33) 


(4) 


4.13 


4.69 


12.5 


100 


253 


32.2 


0.79 


-193 (-146,-46) 


-133 (-111,-22) 


(5) 


4.10 


4.25 






253 


30.2 


0.79 


-193 (-150,-43) 


-120 (-101,-19) 


(6) 


4.01 


5.48 






264 


33.6 


0.76 


-226 (-170, -56) 


-153 (-129,-24) 



Table 4: Neutron-star radius, Rmax (iii km), the central density, rif. (in fm and the ratio 
of the maximum neutron-star mass to the solar mass, Mmax/MQ in SU(6) symmetry. 





Rmax 


ric 


Mmax/MQ 


(1) 


11.5 


1.03 


1.57 


(2) 


10.7 


1.29 


1.36 


(3) 


11.4 


1.08 


1.50 


(4) 


11.3 


1.09 


1.52 


(5) 


11.2 


1.10 


1.62 


(6) 


11.8 


0.97 


1.67 


(7) 


11.9 


0.94 


1.52 
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Table 5: Coupling constants, /7rB/v47r, QmbN^ t^ and kmr = fMn /QMB. 
parameters, A^, in the ESC08 model (taken from lRijken et al.l (120 lOf )). 



and the cutoff 



M 


MNN 


MAA 


MSE 




AB(MeV)'^ 


IT 


f 


0.2675 





0.1903 


-0.0772 


1312 


UJ 


9 


a 


2.8158 


2.8158 


2.0863 


1576 




K 


-0.2276^ 


-1.2030 


-0.0926 


-2.1111 




P 


9 


0.6918 





1.3837 


0.6918 


969 




K 


5.6805 





2.4656 


-0.7489 





"The Lu-N vector coupling constant is determined so as to satisfy the nuclear saturation condition (see Table[6]). 

The tensor coupling is then fixed using k,u:n = —0.2276, which is the value given in the ESC08 model. 
''The form factor in the ESC08 model is given by a gaussian function. The present values are therefore the 

cutoff parameters (given in the ESC08 model) multiplied by ^/2, because we use the dipole function, Eq.([71). 

We assume that the cutoff parameter at the a-BB vertex is identical to that at the tt-BB vertex. 



Table 6: Coupling constants (for the ESC08 model) and calculated properties of symmetric 
nuclear matter at uq. The values of K^j, a^^, and are in MeV, and, in the parenthesis 
of the self-energies, the left (right) number is for the Hartree (Fock) contribution. Note that 
the space component of the vector self-energy is again very small as in the case of SU(6) 
symmetry. 



case 


^ctAT 

A-K 


9tN 
Att 


92 (fm-i) 


^3 




04 


Mjv 




ys 




yO 


(a) 


6.02 


6.13 


22.6 


8.51 


249 


20.0 


0.78 


-204 


(-204, 0) 


-145 


(-145,0) 


(b) 


3.22 


5.16 


13.9 


62.5 


249 


36.2 


0.78 


-204 ( 


-120,-85) 


-144 ( 


-122,-22) 


(c) 


3.12 


4.92 






249 


35.3 


0.78 


-204 ( 


-120,-84) 


-135 ( 


-116,-19) 


(d) 


3.27 


5.86 






261 


35.5 


0.75 


-233 ( 


-143,-90) 


-163 ( 


-138,-25) 



Table 7: The a-Y couphng constants. For details, see the text. 







(A) 






(B) 




case 




9lJ^^ 


9lJ^^ 


9^/4^-^ 


9l^/4rr 


9ls/4^ 


(a) 


5.32 


2.63 


2.71 


2.92 


1.05 


1.45 


(b) 


7.17 


3.32 


3.63 


3.47 


1.02 


1.70 


(c) 


6.59 


3.01 


3.32 


3.09 


0.85 


1.50 


(d) 


5.22 


3.04 


2.50 


2.77 


1.14 


1.29 



-36- 



Table 8: Neutron-star radius, Rmax (in km), the central density, ric (in fm^'^), and the ratio 
of the maximum neutron-star mass to the solar mass, M^ax/^Q, in the ESC08 model. For 
details, see the text. 







(A) 




(b; 




case 


^max 


ric M^ax/MQ 


■^max 


ric 


Mmax/MQ 


(a) 


11.2 


1.04 2.00 


11.6 


0.98 


1.84 


(b) 


11.2 


1.09 1.83 


11.5 


1.05 


1.62 


(c) 


11.2 


1.06 1.86 


11.6 


1.02 


1.71 


(d) 


11.5 


1.00 1.93 


11.7 


0.99 


1.77 



